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Abstract 

Accurate rate coefficients for molecular vibrational transitions due to collisions with H 2 , critical 
for interpreting infrared astronomical observations, are lacking for most molecules. Quantum cal¬ 
culations are the primary source of such data, but reliable values that consider all internal degrees 
of freedom of the collision complex have only been reported for H 2 -H 2 due to the difficulty of the 
computations. Here we present essentially exact full-dimensional dynamics computations for rovi- 
brational quenching of CO due to H 2 impact. Using a high-level six-dimensional potential surface, 
time-independent scattering calculations, within a full angular-momentum-coupling formulation, 
were performed for the deexcitation of vibrationally excited CO. Agreement with experimentally- 
determined results confirms the accuracy of the potential and scattering computations, representing 
the largest of such calculations performed to date. This investigation advances computational quan¬ 
tum dynamics studies representing initial steps toward obtaining CO-H 2 rovibrational quenching 
data needed for astrophysical modeling. 
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Quantum mechanical studies of inelastic processes in molecular collisions began with the 
development of the nearly-exact, close-coupling (CC) method for rotational transitions in 
atom-diatom collisions by Arthurs and Dalgarno [T] and Takayanagi [2] in the 1960s. Over 
the past hve decades tremendous advances in computational processing power and numerical 
algorithms have allowed high-level computations of inelastic processes as well as reactive 
collisional dynamics of large molecular systems, the latter using primarily time-dependent 
approaches |5H7|. Until recently, however, the largest full-dimensional inelastic studies for 
a system of two colliding molecules have been limited to H 2 -H 2 collisions in six-dimensions 
(6D), which were performed with both time-independent [S] and time-dependent [HI [ID] 
approaches for solving the Schrodinger equation. However, to make these computations 
possible, the authors resorted to various angular-momentum decoupling approximations with 
uncertain reliability. It was only recently [IIHID] that these decoupling approximations 
were relaxed and essentially exact full-dimensional CC computations for the H 2 -H 2 system 
became feasible. Despite the large internal energy spacing of H 2 , which allows the basis 
sets to be relatively compact, these calculations are computationally demanding. Replacing 
one H 2 molecule with a molecule that has smaller internal energy spacing, such as carbon 
monoxide, would further increase the computational demands. Whether the CC method 
could in practice be used to describe such a diatom-diatom system in full-dimensionality 
remains an open question. 

First detected in the interstellar medium in 1970 [H], carbon monoxide is the second 
most abundant molecule, after H 2 , in most astrophysical environments. CO has been the fo¬ 
cus of countless theoretical astrophysical studies and observations, being detected in objects 
as distant as high redshift quasars [15] to cometary comae in our solar system [16] to the 
atmospheres of extrasolar giant planets m- Most studies have focused on pure rotational 
transitions observed in the far infrared to the radio or electronic absorption in the near 
ultraviolet. Over the past decade, however, near infrared (NIR) emission of CO due to the 
fundamental vibrational band near 4.7 /rm has been detected in a variety of sources including 
star-forming regions in Orion with the Infrared Space Observatory [TB] and protoplanetary 
disks of young stellar objects [I9H2I] with the Gemini Observatory and the Very Large Tele¬ 
scope. In addition, pure rotational transitions, but in the hrst two vibrationally excited 
states (ui = 1 and 2), were detected by the Submillimeter Array in the circumstellar shell of 
the much-studied evolved star IRC-l-10216 [22]. In particular, high resolution observations 
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of the CO fundamental band probe the physical conditions in the inner disk region, ~10-20 
astronomical units (AU), the site of planet formation. Detailed modeling of such environ¬ 
ments requires state-to-state inelastic rovibrational excitation rate coefficients for CO due 
to H 2 collisions, but current simulations are limited to approximate scaling methods due to 
the dearth of explicit data [IHl 1^ - 

In this Article, we address the two issues outlined above by advancing the state-of-the- 
art for inelastic quantum dynamics with a full-dimensional investigation using an accurate 
potential energy surface relevant to this scattering process, with a particular emphasis on the 
important region of the van der Waals complex. This was made possible through the accurate 
computation and precise htting of a 6D CO-H 2 potential energy surface (PES) in the relevant 
region of the formaldehyde tetra-atomic system and the further development of the inelastic 
diatom-diatom time-independent scattering code, TwoBC [21], which performs full angular- 
momentum coupling, the CC formalism [1], including vibrational degrees-of-freedom. We 
hrst briefly describe the new CO-H 2 PES and its testing through comparison of rotational 
excitation calculations using the 6D PES and 4D PESs (neglecting vibrational motion) and 
available experimental results. The full-dimensional (6D), essentially exact, computations 
for rovibrational quenching of CO(ui = 1) due to H 2 collisions are presented resolving a 
two-orders-of-magnitude discrepancy between earlier 4D calculations which adopted various 
approximations [251128] . Finally, the current results are consistent with the rovibrational 
quenching measurements for the CO-H 2 system, performed at the Oxford Physical Chemistry 
Laboratory from 1976-1993, which no prior calculation has yet been able to adequately 
explain [29H3T]. 


RESULTS 

The CO-H2 potential energy snrface. The CO-H 2 interaction has been of consider¬ 
able interest to the chemical physics community for many decades with one of the hrst 4D 
surfaces for the electronic ground state constructed by Schinke et al. [32], which was later 
extended by Bacic et al. [25 . The group of Jankowski, Szalewicz, and coworkers performed 
accurate 5D and 6D electronic energy calculations, but averaged over monomer vibrational 
modes, and performed several hts to obtain a series of 4D rigid-rotor surfaces, referred to 
as the V98 [33], V04 [M], and V12 [35] PESs. A 6D PES for formaldehyde was constructed 
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earlier by Zhang et al. [36], but it was developed for reactive scattering applications and 
consequently, limited attention was given to the long-range CO-H 2 van der Waals conhg- 
uration. Therefore as a prerequisite to 6D inelastic dynamics studies, we carried out an 
unprecedented potential energy calculation including over 459,756 energy points (see Meth¬ 
ods for details). The potential energy data were then £t using the invariant polynomial 
method with Morse-type variables in terms of bond-distances [371138] , The resulting 6D 
PES, referred to as V6D, is shown in Fig. [TJfor one sample conhguration. Some features of 
V6D are also illustrated in Supplementary Figs. S2 - S4. 

Cross sections and rate coefficients. Time-independent quantum scattering calcula¬ 
tions were performed using the CC formulation of Arthurs and Dalgarno [1] as implemented 
for diatom-diatom collisions in the 4D rigid-rotor approximation in MOLSCAT [3^ and ex¬ 
tended to full-dimensional dynamics as described by Quemener, Balakrishnan, and coworkers 
[nidsi in TwoBC. In the first set of scattering calculations, the new 6D PES was tested for 
pure rotational excitation from CO(ni = 0, ji = 0,1), where Vi and ji are the vibrational 
and rotational quantum numbers, respectively. The crossed molecular beam experiment of 
Antonova et al. m, who obtained relative state-to-state rotational inelastic cross sections, 
is used as a benchmark. The experimental cross sections were determined at three center- 
of-mass kinetic energies (795, 860, and 991 cm“^), but with an initial state distribution of 
CO estimated to be 75 ± 5% for ji = 0 and 25 ± 5% for ji = 1. Antonova et al. normalized 
the relative cross sections to rigid-rotor calculations done with MOLSCAT using the V04 
PES of Jankowski and Szalewicz [33] . Comparison of the experiment with new 4D rigid- 
rotor calculations on V12 and full-dimensional calculations on V6D are shown in Fig. We 
find no difference in the excitation cross sections when using V04 or V12, while the RMS 
cross section errors between the normalized experiment and the V12 and V6D calculations 
range from 0.56-0.89x10“^® cm^ and 0.55-0.95x10“^® cm^, respectively (See Supplementary 
Table SI and Supplementary Note 1). Clearly, a 4D rigid-rotor treatment of the dynamics 
is sufficient for describing rotational excitation at these relatively high energies. 

The importance of full dimensionality for rotational excitation becomes more evident 
as the collision energy is reduced (see also Supplementary Figs. S5 and S6). Low-energy 
excitation cross sections for the process 

CO(ni = 0, ji = 0) -|- 112(^2 = 0,^2 = 0) —)■ CO(n( = 0, = 1) -|- 112(^2 = 0) J 2 = 0); (1) 
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or (0000)—j-(0100), using the notation defined in Methods, were measured by Chefdeville 
et al. [H] in a crossed-beam experiment. They obtained the excitation cross section for 
center of mass kinetic energies from 3.3 to 22.5 cm“^. Though their energy resolution was 
limited, three broad features were detected near 6, 13, and 16 cm“^ attributed to orbiting 
resonances. Computed cross sections for process ([^ using the 4D V12 PES and the full¬ 
dimensional V6D PES are presented in Fig. [^. While both computations reveal numerous 
resonances, the resonances are generally shifted by 2-3 cm“^ between calculations. The 
energy and magnitude of the resonances are very sensitive to the details of the PESs, but 
differences may also be due to the relaxation of the rigid-rotor approximation with the use 
of V6D in TwoBC. In Fig. 3b, the two calculations are convolved over the experimental 
energy resolution and compared to the measured relative cross sections. Except for the 
peak near 8 cm“^, the current 6D calculation appears to reproduce the main features of the 
experiment. RMS errors are found to be 0.355 and 0.228x10“^® cm^ for the V12 and V6D 
PESs, respectively. Further details on the rotational excitation calculations can be found in 
Supplementary Note 1. 

Now that the improved performance of the V6D potential for pure rotational excitation 
is apparent, we turn to rovibrational transitions. As far as we are aware, prior experimental 
[2nHSIl 1121 US] and theoretical [2SH2H] studies are limited to the total quenching of CO(ni = 
1). In the scattering calculations of both Bacic et al. [2S1 [2S] and Reid et al. [27], the 4D 
potential of Bacic et al. ESI was adopted, in which two coordinates were fixed (r 2 = 1.4 ao, 
0 = 0), and various combinations of angular-momentum decoupling approximations for 
the dynamics were utilized (e.g., the infinite order sudden, lOS, and coupled-states, CS, 
approximations; see Supplementary Note 2). More recently. Flower [28] performed CC 
calculations on a parameterization of the 4D PES of Kobaysashi et al. jHj. These four 
sets of 4D calculations for quenching due to para-H 2 (j 2 = 0) are compared in Fig. |^to the 
current 6D/CC calculations on the V6D surface for the case of j 2 = J 2 = 0- Corresponding 
state-to-state and total cross sections for collisions with ortho-H 2 are given in Supplementary 
Figs. S7-S9 with further details in Supplementary Note 2 

A cursory glance at Fig. reveals a more than two orders of magnitude discrepancy 
among the various calculations. The large dispersion for the previous calculations is due to 
a combination of reduced dimensionality and decoupled angular momentum which makes it 
difficult to assess the reliability of each approximation. The current results, however, remove 
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these uncertainties by utilizing i) a full-dimensional (6D) PES, ii) full-dimensional (6D) 
dynamics, and iii) full angular-momentum coupling. The sharp peaks in the cross sections 
over the 1-10 cm“^ range in the 6D/CC results are due to resonances |1S1 SH] supported by 
the CO-H 2 van der Waals potential well. These resonances, in the rovibrational quenching 
of CO by H 2 , are reported here for the hrst time and their prediction is made possible 
through our high-level treatment of the dynamics. Computations for transitions involving 
other initially excited ^2 states and inelastic H 2 (^2 7^ transitions are presented and 
compared in Supplementary Figs. SIO and Sll and Supplementary Note 2. 

Using the current 6D/CC cross sections and the 4D/IOS-CS results of Reid et al, rate 
coefficients for a Maxwellian velocity distribution are computed and compared in Fig. |^to 
total deexcitation measurements [29H3T] reported by Reid et al. m The comparison is not 
straight-forward because i) the measurements correspond to an initial thermal population 
of H 2 rotational states, ii) the initial rotational population of CO(ui = 1, ji) was unknown, 
iii) the experimental rate coefficients for ortho-H 2 are estimated from para- and normal-H 2 
measurements, and iv) the contribution from a quasi-resonant channel, 

CO(ui = 1) + H2 (u 2 = 0, j2 = 2) ^ COK = 0) + H2 (u' = 0, j' = 6) + AF, (2) 

dominates the para-H 2 case for T > 50 K. Fig. displays the CO(ui = 1) rovibrational 
deexcitation rate coefficients from the current 6D/CC calculations for collisions of ortho-H 2 
for ^2=1 and 3, separately. These rate coefficients are summed over all j[ and j 2 =^, 3, and 
5. Assuming a Boltzmann average at the kinetic temperature T of the H 2 rotational levels 
j 2 = 1 and 3, presumed to be representative of the experimental conditions, rate coefficients 
are computed and found to be in good agreement with the measurements above 200 K. From 
Fig. [^, the contribution from ^2=3 is seen to be only important above 150 K. The remaining 
difference with experiment at low temperatures may be due to the fact mentioned above 
that the ortho-H 2 rate coefficients are not directly measured, but deduced from normal-H 2 
and para-H 2 experiments. In particular, Reid et al. assume that the ortho/para ratio in the 
normal-H 2 measurements is 3:1, i.e. statistical, and independent of temperature. Further, as 
stated above, the experimental CO rotational population distribution in ui = 1 is unknown. 
Nevertheless, the current 6D/CC computations are a signihcant advance over the 4D results 
of Reid et al. which also correspond to a Boltzmann average of rate coefficients for j 2 =l 
and 3. 
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As indicated above, the situation for para-H 2 collisions is more complicated due to the 
quasi-resonant contribution (|^, a mechanism not important for ortho-H 2 . Boltzmann- 
averaged rate coefficients are presented in Fig. including ^2=0 and 2 summed over ^ 2 = 
0, 2, and 4, with and without the quasi-resonant contribution, j 2 = 2 —)■ = 6. While 
the current 6D/CC results and the 4D calculations of Reid et al. are in agreement that 
the quasi-resonant contribution becomes important for T > 50 K, the relative magnitude 
compared to the non-resonant transitions from the 6D/CC calculation is somewhat less 
than obtained previously with the 4D potential. This is partly related to the fact that the 
6D/CC rate coefficients for ^2=0 are signihcantly larger than those of Reid et al. (see also 
the corresponding cross sections in Fig. |^and in the Supplementary Fig. S9). Compared to 
the experiment, we obtain excellent agreement for T ^ 150 K, but are somewhat smaller at 
higher temperatures. This small discrepancy may be related to the unknown 00(^1 = 1, ji) 
rotational population in the measurement. Nevertheless, it is only the 6D/CC computa¬ 
tions, i.e., dynamics in full dimensionality with full angular momentum coupling, which can 
reproduce the measurements for both ortho- and para-H 2 . Further, the computations for the 
quasi-resonant process (2) were the most challenging reported here due to the requirement 
of a very large basis set (see Supplementary Table S2) which resulted in long computation 
times, a large number of channels, and usage of signihcant disk space (~0.5 TB per partial 
wave). In total, the cross sections given here consumed >40,000 CPU hours. 


DISCUSSION 

The current investigation of the CO-H 2 inelastic collision system has been performed 
with the intent of minimal approximation through the computation of a high-level poten¬ 
tial energy surface, robust surface htting, and full-dimensional inelastic dynamics with full 
angular-momentum coupling. That is, within this paradigm for studying inelastic dynamics, 
we have advanced the state-of-the-art for diatom-diatom collisions through this unprece¬ 
dented series of computations. The approach has been benchmarked against experiment for 
pure rotational and rovibrational transitions giving the most accurate results to date within 
the experimental uncertainties and unknowns. The accuracy and long-range behavior of the 
6D potential energy surface is found to be comparable to previous, lower-dimensional sur¬ 
faces. The agreement of the current computation for the CO(ni = 1) rovibrational quenching 
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with measurement resolves a long-standing (more than two decades) discrepancy and jus- 
tihes the requirement of a full-dimensional approach. This methodology can now, though 
with signihcant computational cost, be applied to a large range of initial rotational levels for 
ui = 1 and for higher vibrational excitation to compute detailed state-to-state cross sections 
unobtainable via experiment. 

This advance in computational inelastic scattering is particularly timely as ground-based 
(e.g., the Very Large Telescope (VLT)) observations have focused on CO rovibrational emis¬ 
sion/absorption in a variety of astrophysical objects, while related observations are in the 
planning stages for future space-based telescopes (e.g., NASA’s James Webb Space Tele¬ 
scope). In particular, we are in an exciting era of investigation into the properties of pro¬ 
toplanetary disks (PPDs) around young stellar objects |1^. PPDs provide the material for 
newly-forming stars and fledgling planets. The CO fundamental band (|Aui| = 1) is a tracer 
of warm gas in the inner regions of PPDs and, with appropriate modeling, gives insight into 
disk-gas kinematics and disk evolution in that zone where planets are expected to be forming. 
A recent survey of 69 PPDs with the VLT [19] detected CO vibrational bands in 77% of the 
sources including ui = 1 —0, Ui = 2 —)■ 1, ui = 3 —)■ 2, and even ui = 4 —)■ 3 in a few cases. 
Remarkably, rotational excitation as high as ji = 32 was observed. However, the modeling 
of PPDs, and other astrophysical sources with CO vibrational excitation, is hindered by the 
lack of rate coefficients due to H 2 collisions. We are now in an excellent position to provide 
full-dimensional, state-to-state CO-H 2 collisional data which will not only have a profound 
impact on models characterizing these intriguing environments that give birth to planets, 
but also aid in critiquing current theories used to describe their evolution. 


METHODS 

Potential energy computations. The potential energy computations were performed 
using the explicitly correlated coupled-cluster (CCSD(T)-F12B) method jlHHH], as imple¬ 
mented in MOLPRO2010.1 (50] • The cc-pcvqz-fl2 orbital basis sets [SI] that have been 
specifically optimized for use with explicitly correlated F12 methods and for core-valence 
correlation effects have been adopted. Density fitting approximations [19] were used in all ex¬ 
plicitly correlated calculations with the AUG-CC-PVQZ/JKFIT and AUG-CC-PWCVQZ- 
MP auxiliary basis sets [52l [53]. The diagonal, fixed amplitude 3G(FIX) ansatz was used. 


which is orbital invariant, size consistent, and free of geminal basis set superposition error 
(BSSE) [nUin^. The default CCSD-F12 correlation factor was employed in all calculations 
and all coupled cluster calculations assume a frozen core (C:ls and 0:1s). The counter-poise 
(CP) correction [56] was employed to reduce BSSE. Even though the explicitly correlated 
calculations recover a large fraction of the correlation energy, the CP correction is still neces¬ 
sary, mainly to reduce the BSSE of the Hartree-Fock contribution. Benchmark calculations 
at the CCSD(T)-F12b/cc-pcvqz-fl2 level were carried out on selected molecular configura¬ 
tions and results were compared with those from the conventional CCSD(T) method using 
aug-cc-pV5Z and aug-cc-pV6Z basis sets. Results showed that the CP corrected interaction 
energy agrees closely with those derived from CCSD(T)/aug-cc-pV6Z. 

To construct the potential energy surface (PES), the computations were performed on a 
six-dimensional (6D) grid using Jacobi coordinates as shown in Supplementary Fig. SI. R 
is the distance between the center-of-mass of CO and H 2 . ri and r 2 are the bond lengths of 
CO and H 2 , respectively. 61 is the angle between ri and R, 62 the angle between r 2 and R, 
and 0 the out-of-plane dihedral or twist angle. In the potential energy computations, the 
bond lengths are taken over the ranges 1.7359 < ri < 2.5359 ao and 1.01 < r 2 < 1.81 ao, 
both with a step-size of 0.1 ao. For R, the grid extends from 4.0 to 18.0 ao with step-size of 
0.5 ao for R < 11.0 ao and 1.0 ao for R > 11.0 ao. All angular coordinates were computed 
with a step-size of 15° with 0 < < 360° and 0 < 02,0 < 180°. Additional points were 

added in the region of the van der Waals minimum. 


The PES fit. The CO-H 2 interaction PES has been fitted in 6D using an invariant 
polynomial method mm- The PES was expanded in the form, 

N 

nvi • • • Be) = E c„.-n,y’;'yTlyTy?vTvf + v?yry?yn (3) 

nl-"n6 

where Ui = is a Morse-type variable. The internuclear distances di between two atoms 

are defined as di = dnu', 02 = dow, 03 = 0ch', 04 = 0ch, 05 = 0oh, and 06 = 0co- The total 
power of the polynomial, N = ni + n 2 + + n 4 + + uq, was restricted to 6. Expansion 

coefficients were obtained using weighted least-squares fitting for potential energies 

up to 10,000 cm“^. The root-mean-square (RMS) error in the fit of the PES is 14.22 cm“^ 
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which included 398,218 different geometries. This RMS error can be compared to that of 277 
cm“^ for the 6D reactive surface of Zhang et al. [SB]. From the computed energy points, the 
global minimum of the total potential corresponds to the collinear arrangement H-H-C-0 
(01 = 0, 02 = 0, 0 = 0) with a depth of -85.937 cm“^ at R = 8.0 ao with ri and r 2 at their 
respectively equilibrium positions. This compares to the global minimum of the interaction 
potential obtained by Jankowski et al. [35] from their fitted V12 PES averaged over ri and 
r 2 : R = 7.911 ao and -93.651 cm“^. Note that this comparison is only suggestive as the 
global minimum in the total and interaction potentials coincide only when bond-lengths (ri 
and r 2 ) are fixed at their equilibrium values as illustrated in Fig. 1 and which is not the 
case for V12. 


Scattering theory and computational details. The quantum scattering theory for 
a collision of an F-state atom with a rigid-rotor was developed la EMS] based on the 
close-coupling (CC) formulation of Arthurs and Dalgarno [1]. Details about its extension to 
diatom-diatom collisions with full vibrational motion can be found in Refs. nans]. In this 
approach, the interaction potential R(R, ri, r 2 , 0i, 02 , 0) is expanded as, 

l/(R,ri,r2,0i,02,0) = Y. ^\i,\2,\i2iri,r2,R)Yx^,x^,xYh,f2,R), (4) 

A1A2A12 

with the bi-spherical harmonic function expressed as, 

Fai,A2,Ai2(d,^2,R) = Y (•^l"iAiA2mA2|Ai2mAi2) 
mx^mx2'fnxi2 

xFAimA^ (5) 

where 0 < Ai < 10, 0 < A 2 < 6 was used in the scattering calculations. Due to the symmetry 
of H 2 , only even values of A 2 contribute. 

For convenience, the combined molecular state (CMS) notation is applied to describe 
a combination of rovibrational states for the two diatoms. A CMS represents a unique 
quantum state of the diatom - diatom system before or after a collision. The CMS will be 
denoted as (nijin2j2)- j are the vibrational and rotational quantum numbers. 

The rovibrational state-to-state cross section as a function of collision energy E is given 
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by, 


'^VijiV2j2->-V 



X 



( 6 ) 




where (^ 171 ^ 2 ^ 2 ) and are, respectively, the initial and hnal CMSs of CO-H 2 , the 

wave vector k'^ = 2fj.E/h‘^, and S is the scattering matrix. I is the orbital angular momentum 


and J the total collision system angular momentum, where J = 1 + ji2 and ji2 = ji + h- 

Thorough convergence testing was performed in the scattering calculations by varying all 
relevant parameters. The CC equations were propagated for each value of R from 4 to 18.0 ao 
using the log-derivative matrix propagation method of Johnson [60] and Manolopoulos [ 6 T] , 
which was found to converge for a radial step-size of AR = 0.05 ao. The convergence tests 
of the Ui = 1 —)■ 0 vibrational quenching cross section of CO with respect to the number 
of = 1 rotational channels found that at least 13-15 channels have to be included in 
the Ui = 1 basis set, especially for low-energy scattering. Based on convergence tests with 
respect to the adopted maximum R for the long range part of the PES, we found that 
the cross sections are converged down to the lowest collision energy of 0.1 cm“^. This 
value also guarantees that the rate coefficients are converged for temperatures greater than 
1 K. The number of discrete variable representation points and N^^', the number of 
points in 9i and 6*2 for Gauss-Legendre quadrature, and and the number of points 
in cj) for Gauss-Hermite quadrature, N^, which were applied to project out the potential 
expansion coefficients were all tested for convergence with the hnal adopted values given in 
Supplementary Table S2. The basis sets and the maximum number of coupled channels are 
also presented in Supplementary Table S2. 

The resulting integral cross sections were thermally averaged over a Maxwellian kinetic 
energy distribution to yield state-to-state rate coefficients as function of temperature T, 



'OO 


^VljlV2j2^v'J[v'^j'^ (T) 


(E) exp(-/3E)dE, (7) 


Ea, 


VljlV2j2^v[j[v!2j!2 


where m is the reduced mass of the GO-H 2 complex, (3 = (ksT) and ks is Boltzmann’s 
constant. 
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FIG. 1. The CO-H 2 interaction potential energy surface V6D. The potential surface is 
constructed in the 6D diatom-diatom Jacobi coordinates (i2, ri, r 2 , , 02 j </>), and r 2 are bond 

lengths, R the internuclear distance between CO and H 2 center of masses, 61 and 62 the angles 
between R and ri and r 2 , and (j) the dihedral or twist angle. See Supplemental Fig. 1. Here the 
dependence of the potential surface on bond lengths rco = ri and thj = r 2 is shown with R = 8 
ao, 01 = 180°, 02= 0, and (j)= 0. Note that the CO(ri) and H 2 (r 2 ) diatom potentials have been 
subtracted. 
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FIG. 2. State-to-state cross sections for rotational excitation of CO(?;i = 0,ji = 0,1) 
by collisions with H 2 . Theoretical results of full-dimensional calculations on V6D and 4D rigid- 
rotor calculations on VI2 are compared with normalized experimental results m for collision 
energies of (a) 795 cm~^, (b) 860 cm“^, and (c) 991 cm“^. The calculations were performed 
for H 2 (u 2 = 0, j 2 = 0), but the experimental H 2 rotational distribution was undetermined. The 
error bars correspond to twice the estimated standard deviation in the weighted means of the 
measurements m- 
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FIG. 3. Low energy excitation cross sections. = 0 —1 cross sections for CO(?;i = 0) due 
to collisions with H 2 (u 2 = 0 ,j 2 = 0) are shown as a function of collision energy, (a): computed 
cross sections using the 4D V12 and 6D V6D PESs. (b): computed cross sections convolved over 
the experimental beam energy spread (lines) compared to the relative experiment of Chefdeville et 
al. |41j (circles with error bars). The error bars on the experimental cross sections of Chefdeville 
et al. represent the statistical uncertainty at a 95% confidence interval. 
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FIG. 4. Total theoretical cross sections for the vibrational deexcitation of CO(z;i = 1) 
by para-H 2 . Current 6D/CC results are compared to previous 4D calculations. The 4D results of 
Bacic et al. [26] and Reid et al. m do not distinguish CO rotational states, while the 4D results 
of Bacic et al. |25|, 4D results of Flower |28|, and the current 6D/CC results are for initial ji = 0 
summed over all final CO(ui = 0, j^). In every case, the H 2 rovibrational state remains unchanged, 
V2,j2 = V2,j2 = 
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FIG. 5. Rate coefficients for the vibrational deexcitation of CO(?;i = 1) dne to H 2 . 

Current 6D/CC (solid, dashed, and dotted lines) and Reid et al. [57] 4D (dot-dashed line) cal¬ 
culations are compared to the total CO(ui = 1) rovibrational quenching experiment (symbols) 
[29I - I3T] . (a) Ortho-H 2 rate coefficients for initial state-resolved, (100j2) —)• (Oj^Oj^), summed over 
j'l and J 2 (dashed and dotted lines) and for a Boltzmann average of initial H 2 (j 2 = 1) 3) (solid and 
dot-dashed lines), (b) Para-H 2 rate coefficients for a Boltzmann average of initial H 2 (j 2 = 0,2) 
with and without the quasi-resonant process Q. Note that the experimental uncertainties are 
smaller than the symbol sizes. 
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Supplementary Information 


SUPPLEMENTARY FIGURES 



FIG. SI. The six-dimensional Jacobi coordinates for CO-H 2 . ri and r 2 are bond lengths, R the 
internuclear distance between CO and H 2 center of masses, 9i and 62 the angles between R and fl 
and r 2 , and (f> the dihedral or twist angle. 
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FIG. S2. Slices of the CO-H 2 interaction potential V 6 D. R dependence of the interaction potential 
V 6 D for representative slices with bond lengths fixed as indicated and ( 6 * 1 , 62 , 4))={0°, 0°, 0°), 
(180°, 0°, 0°), and (90°, 90°, 90°). V 6 D fit (lines), computed ab initio energy points (symbols). 
This figure may be compared to Fig. 2 of Ref. [ST] . 
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FIG. S3. Percent differences between CO-H 2 PESs. Percent difference between the fitted PES and 
ab initio energy points of the interaction potential V6D for representative slices with bond lengths 
fixed as indicated and (0i, 62 , (p)={0°, 0°, 0°), (180°, 0°, 0°), and (90°, 90°, 90°). The middle 
panel also displays the percent difference between the V12 PES and the current ab initio data. 
The percent differences only become large in the asymptotic limit as the PES itself approaches 
zero. Also, displayed in the middle panel is the difference between the V12 PES and the current 
ab initio data. It is seen that there is good agreement, even in the long-range limit, except for 
the 9i = 180° configuration. This difference may partly be due to the fact that V12 is a 4D PES 
obtained by averaging a 6D PES over the CO and H 2 ground state vibrational wave functions. 
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FIG. S4. Contour plot of the CO-H 2 V6D PES. The contour plot is given as a function of Q\ and 
02- .R = 7-5 ao, (j) = 0°, ri = 2.1359 ao, r 2 = 1.4011 uq. The legend gives the contour energies in 
cm“^. This figure may be compared to Fig. 1 of Ref. |Slj . 
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FIG. S5. CO rotational excitation cross sections due to H 2 collisions. Rigid-rotor approximation 
= 0 —)• j'l excitation cross sections of CO(r;i = 0,ji) by collisions with H 2 (r ;2 = 0,j2 = 0) 
using the V6D (EE and VV) and V12 potentials compared to the full-dimensional calculation. See 
Supplementary Note 1 for further details. Collision energy: (a) 50 cm“^, (b) 100 cm“^, and (c) 
500 cm“^. 
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FIG. S6. CO rotational excitation cross section difference for H 2 collisions. Difference of ji = 
0 —)• j'l excitation cross sections of CO(ni = 0,ji) by collisions with H 2 (n 2 = 0,^2 = 0) using the 
V6D (EE and VV), V12 potentials and those of the full-dimensional calculation (Full-D, V6D). 
See Supplementary Note 1 for further details. Collision energy: (a) 50 cm~^, (b) 100 cm~^, and 
(c) 500 cm“^. 
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FIG. S7. C0-para-H2 vibrational quenching state-to-state cross sections. State-to-state cross 
sections for the vibrational quenching of CO from (ui = 1, ji = 0) to = 0, jj), j[=0, 2, 4, • • •, 
24, due to para-H 2 (u 2 = 0,^2 = 0) collisions, or in CMS notation (1000)—)- (OjJOO). 
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FIG. S8. CO-ortho-H 2 vibrational quenching state-to-state cross sections. State-to-state cross 
sections for the vibrational quenching of CO from (m = 1, ji = 0) to (?;( = 0, j(), j[=0, 2, 4, • • •, 
24, by ortho-H 2 (n 2 = 0,j2 = 1), or in CMS notation (1001)—)- (Oj(Ol). 
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FIG. S9. Total CO-ortho-H 2 vibrational quenching cross sections. Total cross sections for the 
vibrational quenching of CO(ui = 1) + H 2 (u 2 = 0, j 2 = 1) to CO(u'^ = 0) + ortho-H 2 (u 2 = 0, J 2 =l). 
Similar to Fig. 4 of the main text, except for ortho-H 2 . 
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FIG. SIO. Total CO-para-H 2 vibrational quenching cross sections for different H 2 states. Total 
cross sections for the vibrational quenching of CO from initial states (1000) (top panel) and (1002) 
(bottom panel) to CO(n( = 0) + para-H2(j2=0, 2, 4, 6). 
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FIG. Sll. Total CO-ortho-H 2 vibrational quenching cross sections for different H 2 states. Total 
cross sections for the vibrational quenching of CO from initial states (1001) (top panel) and (1003) 
(bottom panel) to CO(n(=0) + ortho-H 2 (j 2 =l, 3, 5). 
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SUPPLEMENTARY TABLES 


TABLE SI. CO-H 2 rotational excitation errors. RMS error (10 cm^) of the three rigid-rotor 
calculations and V6D in comparison with experiment |S2] at E = 795, 860, and 991 cm“^. 


E (cm 

V6D-full 

V6D-EE 

V6D-VV 

V12 

795 

0.549 

0.546 

0.552 

0.556 

860 

0.571 

0.598 

0.572 

0.551 

991 

0.946 

0.989 

0.946 

0.892 


12 





TABLE S2. Parameters used in the CO-H 2 scattering calculations. 


Basis se1‘ 

1 

Ne, 

Ne2 


Nr, 

Nr, 

Ai 

A 2 

No. of channels 

4D rotational calculation 

TwoBC 

ji = 30, 32 

= 2 

14 

14 

8 

1 

1 

8 

4 

1344 

MOLSCAT 

(M 

o' 

CO 

II 

= 2 

14 

14 

8 



8 

4 

1344 

6D TwoBC calculation 

para-H 2 -CO 

[(0,30;1,20)(0,2)] 

14 

14 

8 

20 

20 

10 

6 

4258 

para-H 2 -CO 

[(0,22;1,20)(0,4)] 

14 

14 

8 

20 

20 

10 

6 

7341 

para-H 2 -CO 

[(0,20;1,15)(0,6)] 

12 

12 

8 

18 

18 

8 

4 

9832 

ortho-H 2 -CO 

[(0,30;1,20)(0,3)] 

14 

14 

8 

20 

20 

10 

6 

5632 

ortho-H 2 -CO 

[(0,22;1,20)(0,5)] 

14 

14 

8 

20 

20 

10 

6 

10251 


^ Basis set [(vi = 0,jvi=o', = 1, jvi=i){v 2 = 0,j.u2=o)] is presented by the maximum rotational quantum 

number jy^ and jy^ included in each relevant vibrational level vi and V 2 for CO and H 2 , respectively. 
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SUPPLEMENTARY DISCUSSION 


Supplementary Note 1: Pure rotational scattering 


The V6D PES was adopted in CC scattering calculations using both TwoBC [S3] and 
MOLSCAT |S4j . Separate calculations were performed with the bond lengths of CO and H 2 
fixed at: 1) equilibrium bond lengths, ri=2.1359 ao and r2=1.4011 ao, denoted as EE; and 
2) vibrationally-averaged bond lengths, ri=2.13992 ao and r 2 =l.448736 ao, denoted as VV. 
The 4-dimensional (4D) V12 potential of Jankowski et al. |S5] . which was constructed by 
averaging over their 6D potential energy data with CO(ni = 0) and H 2 (n 2 = 0) vibrational 
wave functions, was used in MOLSCAT calculations only. For the rigid-rotor scattering 
calculations, TwoBC and MOLSCAT give almost identical CO state-to-state excitation cross 
section using the 4D version of V6D. As an example. Supplementary Fig. displays the 
CO state-to-state excitation cross sections from jJ = 0 for collision energies of 50, 100, and 
500 cm“^ due to H 2 (n 2 = 0,^2 = 0) using the V12 potential with MOLSCAT and the two 


4D (EE and VV) versions of V6D with TwoBC. Supplementary Fig. S6 gives the differences 


in these cross sections compared to V6D (full-dimensional) results. There is generally good 
agreement among all approaches with a small difference depending on which CO and H 2 
bond-lengths are adopted. Comparison is also made to the full 6D calculation which agrees 
best with the rigid-rotor VV results. This conhrms both the validity of the rigid-rotor 
approximation and the suggestion of Jankowski and Szalewicz [S6] that PESs constructed 
with vibrationally-averaged bond lengths are preferred for scattering calculations. Further, 
the discrepancy in the rigid-rotor calculations and the 6D results are seen to increase, at 
least for low j [, with decreasing collision energy supporting the claim in the main text that 
accurate prediction of low-energy dynamics requires full-dimensionality. 

Antonova et al. |S2] measured the relative state-to-state rotationally inelastic cross sec¬ 
tions for excitation of CO by H 2 in a crossed molecular beam experiment at collision energies 
of 795, 860, and 991 cm“^ as discussed in the main text. Comparison between the current 
calculated cross sections and the measurements give a root mean square (RMS) error com¬ 
puted by. 


RMS = 


?(o'expt ^theo)" 
^ N 


(SI) 
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where a is the cross section and N the number of data points. The RMS errors for the 4D 
version of V6D with EE and VV bond lengths, the 4D V12 PES, and the V6D PES are 
shown in Supplementary Table 1 with the cross sections for the latter two shown in Fig. 2 
of the main text. The current V6D results for the final CO j[ distributions (Fig. 1) are seen 
to be in excellent agreement with experiment, validating the accuracy of the htted potential 
as well as the present scattering calculations. 

Chefdeville et al. |S7] reported crossed-beam experiments and quantum-mechanical cal¬ 
culations performed for the CO(ji=0) -|- H2(j2=0) —)■ CO(j(=l) -|- H2(j2=0) process. The 
4D V04 PES of Jankowski and Szalewicz |S6j was used in their scattering calculations. 
To compare with the measurements, the computed cross sections were convolved with the 
experimental kinetic energy distribution according to. 


a(B) = 


1 


Jo 


a{E') exp 


{E-E 


r\2 


2J2 


dE\ 


where 5 is given by 


5 = ^2(0.077 + 0.051E - 3.6 x 


(S2) 


(S3) 


The detailed cross sections and discussion are presented in the main text. 


Supplementary Note 2: Full-dimensional Scattering calculations 


Quenching cross sections 


We have performed full-dimensional calculations of the state-to-state cross sections for 
initial CMSs (100j2), for para-H 2 (j2=0 and 2), and ortho-H 2 (j 2 =l and 3), for collision 


energies ranging from 0.1 to 1000 cm“^. As examples, in Supplementary Figs. ^ and ^ the 
state-to-state quenching cross sections for hnal states (0j(0j2), Ji=0, 2, 4, • • •, 24 are shown 
for para- and ortho-H 2 , respectively, where j 2 = i-®-) elastic in H 2 . The cross sections 

have similar behaviour and display resonances in the intermediate energy region primarily 
due to quasibound states supported by the van der Waals well of the interaction potential of 
the initial channel. Resonances and oscillatory structures, which will be analyzed further in 
a future publication, appear up to about 100 cm“\ roughly the depth of the van der Waals 
minimum. Small |Aji| = \j[ — ji\ transitions dominate the quenching; the cross sections 
generally decreasing with increasing j[ with that for j[ = 24 being the smallest. 
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The current 6 D/CC state-to-state cross sections are summed over ][ (for v[ = 0) to 


compare to the total deexcitation 4D cross sections given by Bacic et al. 


j, Reid 


et al. m, and Flower [Sll] as shown for para-H 2 in Fig. 4 of the main text and here 


for ortho-H 2 in Supplementary Fig. ^ As discussed in the main text there is a two order 
of magnitude dispersion between the various calculations which may partly be explained 
by the differing treatments of angular momentum coupling. In the infinite order sudden 
(lOS) approximation, the internal rotation is neglected (effectively setting j = 0), while the 
coupled-states (CS) approach, which is a centrifugal sudden approximation, allows for the 
dynamical treatment of rotational transitions. Bacic et al. |S 8 j performed both CS and lOS 
calculations, but with j 2 = 0 (i.e., only para-H 2 ), while Bacic et al. [S9] and Reid et al. 
[Sinj used a mixed lOS-CS formalism in which the rotational motion of CO was neglected, 
but that of H 2 retained, allowing for signihcant reduction of basis sets. The cross sections 
from Flower isni are likely too large due to an insufficient CO basis, an effect studied by 
Bacic et al. |S 8 j . In fact, using Flower’s basis set in CC calculations on V 6 D, we obtained 
signihcantly larger cross sections than our converged results, verifying the primary cause for 
the discrepancy. However, in most cases differences are related to the adopted PESs and 
dimensionality of the scattering approaches. 


Quenching rate coefficients 

In the calculations of Bacic et al. [S8l |S9] and Reid et al. EE], the total quenching cross 
section from state (ui = 1 ,^ 2 ) was obtained by summing the state-to-state quenching cross 
sections over the final rotational state of H 2 , 

= 1,42 ^ v[ = 0){E) = ^a(ni = 1,42 ^ v[ = 0,42)(E), (S4) 

32 ' 

allowing for H 2 rotational inelastic transitions (with V 2 = v '2 = and no information on 41 
or i[ as the lOS approximation was adopted for CO in both previous calculations.) By ther¬ 
mally averaging the cross sections over the collision energy, the quenching rate coefficients 
from vi = 1,42 were obtained. 

The total thermal rate coefficients for CO ui = 1 —)■ = 0 vibrational quenching were 
obtained by summing k{yi = 1,42 —)■ = 0)(T) over all initial H 2 rotational states 42 
weighted by their populations assuming a Boltzmann distribution at T. As observed by 
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Millikan and Osbnrg |S12] . para-H 2 is a more efficient collision partner than ortho-H 2 for 
vibrational quenching of CO due to the quasi-resonant transition: 

CO(ni = 1) -|- H 2 (n 2 = 0,j2 = 2) —)• CO(n^ = 0) -|- 112(^2 = OO 2 = 6). Guided by 
the measurements, this quasiresonant process was also theoretically investigated by Bacic 
et al. [S9] and Reid et al. |Sinj requiring cross sections for rotationally excited H 2 with 
inelastic H 2 transitions. Supplementary Figs. |S10| and |S11| compare the ^2 state-to-state 
cross sections from the earlier work with the current results. As shown in the top panel 


of Supplementary Fig. SIO, the current 6D/CC results are in reasonable agreement with 
the 4D/CS-IOS calculations of Bade et al. for the 1000—)-0jj00 and 1000—j-0j(02 cases, but 
both results are nearly an order of magnitude larger than obtained by Reid et al. In the 


bottom panel of Supplementary Fig. |S10| , a large scatter in the results is evident for initial 
state (1002), but generally the cross sections increase with increasing Aj 2 = — J 2 - The 

important exception is that contrary to both the current calculations and the 4D/CS-IOS 
results for Reid et al, Bacic et al. hnd much smaller cross sections for the quasi-resonant 
process. The V6D/CC results shown in that hgure for the quasi-resonant transition were the 
most computationally time-intensive due to the larger size of the H 2 basis and may represent 
a lower limit for the cross section at the highest energies. 


Supplementary Fig. |Sll| gives a similar comparison for ortho-H 2 collisions though a quasi¬ 
resonant process is not operable due to much larger asymptotic energy differences. Other¬ 
wise, the trends are very similar to those noted for para-H 2 collisions, though the current 
V6D/CC results fall somewhat between the earlier two 4D calculations. Note that we did 
not consider the H 2 inelastic channel ^2 = 1 5 as Bade et al. |S9j hnd the cross sections 

to be more than an order of magnitude smaller than the j 2 = 1 3 transition. This is 

in contrast to Reid et al. who hnd the j 2 = 1 —t 5 transition to dominate the quenching 
of (1001). As a consequence, the thermally-averaged rate coefficients for the quenching by 
ortho-H 2 obtained by Reid et al, shown in Fig. 5a of the main text, are signihcantly larger 
than the V6D/CC computation and the experiment. These rate coefficients also include 
quenching from the initial state (1003), for which the cross sections are shown in the bottom 
panel of Supplementary Fig. Sll[ though Reid et al. do not tabulate them. Because of the 
complexity of the system and diherences in the PESs and scattering treatments, it is difficult 
to draw conclusions from one-to-one comparisons of the various calculations. However, the 
main dehciency in the lOS approach is that it does not resolve the hnal CO(j() channels 
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which can play a significant role in the cross sections as illustrated in Supplementary Figs.[S7 
and 
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